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Abstract 

We study the role of fundamental constants in an updated recom- 
bination scenario. We focus on the time variation of the fine structure 
constant a, and the electron mass rrif, in the early Universe, and put 
bounds on these quantities by using data from CMB including WMAP 
5-yr release and the 2dFGRS power spectrum. We analyze how the 
constraints are modified when changing the recombination scenario. 

Time variation of fundamental constants is a prediction of theories that 
attempt to unify the four interactions in nature. Many efforts have been 
made to put observational and experimental constraints on such variations. 
Primordial light elements abundances produced at Big Bang nucleosynthesis 
(BBN) and cosmic microwave background radiation (CMB) are the most 
powerful tools to study the early universe and in particular, to put bounds 
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on possible variations of the fundamental constants between those early times 
and the present. 

Previous analysis of CMB data (earlier than the WMAP five-year release) 
including a possible variation of a have been performed by refs. [H [2|, |3l H] 
and including a possible variation of mg have been performed by refs. O U, |5] . 
In March 2008, WMAP team released data collected during the last five years 
[0]. Moreover, new processes relevant during the recombination epoch have 
been taken into account. Indeed, in the last years, helium recombination 
has been studied in great detail [3 [HI [H [10], revealing the importance of 
these considerations on the calculation of the recombination history. Switzer 
& Hirata [8] presented a multi-level calculation for neutral helium recombi- 
nation including, among other processes, the continuum opacity from H I 
photoionization. They found that at 2; < 2200 the increasing H I abundance 
begins to absorb photons out of the He I 2^p — > l^s line, rapidly accelerating 
He I recombination, which finishes at z ~ 1800. Kholupenko et al [11] have 
considered the effect of the neutral hydrogen on helium recombination and 
proposed an approximated formula to take this effect into account. This has 
enabled its implementation on numerical codes such as Recfast, since the 
complete calculations done by Switzer & Hirata take a large amount of com- 
putational time. Another improvement in the recombination scenario is the 
inclusion of the semi-forbidden transition 2^p l^s, the feedback from spec- 
tral distorsions between 2^p l^s and 2^p — > l^s lines, and the radiative 
line transfer. 

The release of new data from WMAP brings the possibility of updating 
the constraints on the time variation of fundamental constants. In this paper 
we study the variation of a and me in the improved recombination scenario. 
It could be argued that me is not a fundamental constant in the same sense as 
a is and that constraints on the Higgs vaccum expectation value (< f >) are 
more relevant than bounds on me- However, in the recombination scenario 
the only consequence of the time variation of < f > is a variation in me. 

The effect of a possible variation of a and/or me in the recombination 
scenario and in the CMB temperature and polarization spectra has been 
analized previously [121 [ISl [H [S]- Here we focus in the effect of the variation 
of a and me on the improved recombination scenario. 

The recombination equations implemented in Recfast in the detailed 
recombination scenario [T3| including the fitting formulae of [11] are: 
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H{z){l + z)^ = [x^x^n^n - /9h(1 - a;p)e-^'^H-/^^-)CH, (1) 

H[z){l + Z)'^^ = (xHeIlXenH«HeI " /^Hell/He " XHeIl)e-"'^Hel.2i.ArM) ^^^^ 

+ fxHeIlXenH«HeI " ^^^/?Hel(/He " XHeIl)e-'''^Hel..3./fcTM \ ^^^^ ^ 
V 5'HoI,lis / 

(2) 



where 

^ _ l + A'HAH^H(l-a^p) 

^" l + i^H(AH + /3HK(l-Xp)' ^ ^ 

^ ^ 1 + i^HclAHc^H(/Hc - XHelQe^'^P-/^^'^ , . 

""^^ 1 + i^Hel(AHe + Munihe " XHeIl)e'^--/'='^M ' ^ ^ 

^"'^^ ^ 1 + neI/3Ll^H(/He - Xn.u)^^^'^^"^^ ' 

The last term in eq. ^ accounts for the recombination through the 
triplets by including the semi-forbidden transition 2^p l^s. As remarked 
in [E] , ttHei is fitted with the same functional form used for the ctHei singlets, 
with different values for the parameters, so the dependences on the funda- 
mental constants are the same, being proportional to a^m^^^^. The two 
photon transition rates Ah and Ahb depend on the fundamental constants 
as Q^me- The photoionization coefficients (3 are calculated as usual from the 
recombination coefficients (with c standing for H, Hel and Hell), so their 
dependencies are known (see for example ||4j). 

The i^H, -^Hei and -R'nei quantities are the cosmological redshifting of the 
H Lya, Hel 2^p-l^s and Hel 2^p-l^s transition line photons, respectively. 
In general, K and the Sobolev escape probability ps are related through the 
following equations (taking Hel as an example): 

Knei = THrf (6) 

fi'HeI,2ip ^HeI,lis^2ip-lisPS 

T^t _ ^'Hel.lis 1 /„x 

-^Hel - 4HeI „ ' \' ) 

fl'HeI,23p ^HeI,lls^23p-lis^'S 
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where ^Hei,2ip-iis and ^Hei,23p-iis are the Einstein A coefficients of the He I 
2^p-l^s and He I 2^p-l^s transitions, respectively. To include the effect of the 
continuum opacity due to H, based on the approximate formula suggested 
by ref. [11], ps is replaced by the new escape probability Pesc = Ps + Pcon,H 
with 

_ 1 

Pcon,H - — T-, (8) 

1 + aHe7 



and 



^7r3/2o-H,ls(z^HeI,2lp)'^HeI,2ip^^D,2ip(l " ^^p) 



ry - gHoI,2lp - ^ — - , . 

I 0„3/2„ f,. _ \,.2 A,, . n ^ \ 



where o"H,is('^Hei,2ip) is the H ionization cross-section at frequency z/Hei,2i 



p 



and Az/D,2ip = '^Hei,2ip y 2fcB7M/'"^HeC^ is the thermal width of the He I 2^p- 
l^s line. The cross-section for photo-ionization from level n is |16j : 



an{Z,hu) = -—^ — {l + nh)-'gu{n,e) (10) 

where gij{n,e) ~ 1 is the Gaunt-Kramers factor, and = Ti/ {171^.00) is the 
Bohr radius, so crH,is(^'Hei,2ip) is proportional to a~^m~^. 

The transition probability rates AHei,2ip-iis and v4Hei,23p-iis can be ex- 
pressed as follows pT] : 



^^-1 = ^4K^^K + ^2|^,)r (11) 



where Uij is the frequency of the transition, and refers to the initial 
(final) state of the atom. First we will analize the dependence of the bra-ket. 
To first order in perturbation theory, all wavefunctions can be approximated 
to the respective wavefunction of hydrogen. Those can be usually expressed 
as exp(— gr/ao) where oq is the Bohr radius and g is a number. It can be 
shown that any integral of the type of Eq. (iTTl) can be solved with a change 
of variable x = r/ao. If the wave functions are properly normalized, the 
dependence on the fundamental constants comes from the operator, namely 
ri + r2. Thus, the dependence of the bra-ket goes as oq. On the other hand, 
ujij is proportional to the difference of energy levels and thus its dependence 
on the fundamental constants is uJij ~ meO^. Consequently, the dependence 
of the transition probabilities of Hel on a and is 
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(12) 
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Figure 1: Ionization history allowing a to vary with time. From left to 
right, the values of ^ are 1.05, 1.00, and 0.95, respectively. The dotted 
lines correspond to the standard recombination scenario, and the solid lines 
correspond to the updated one. 

In Fig. [1] we show how a variation in the value of a at recombination 
affects the ionization history, moving the redshift at which recombination 
occurs to earlier times for larger values of a. The difference between the 
functions when the two different recombination scenarios are considered, for 
a given value of a, is smaller than the difference that arise when varying the 
value of a. Something similar happens when varying mg. 

With regards to the fitting parameters and fone, since detailed cal- 
culation of their values are not available yet, it is not possible to determine 
the effect that a variation of a or me would have on these new parameters. 
Wong et al [15] have shown that they must be known at the 1% level for 
future Planck data. In this letter, however, we are dealing with the 5 yr 
data from WMAP satellite and this accuracy is not required. To come to 
this conclusion, we have calculated the temperature, polarization and cross 
correlation CMB spectra, allowing the parameters and fene to vary at the 
50% level. We found that for the temperature and polarization spectra, the 
variation is always lower than the observational error (1% for temperature 
and almost 40 % in polarization). The largest variations occur in the cross 
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correlation CMB spectra (Cj^). In this case, we have calculated the ob- 
servational errors divided by the value of the Cis of all measured Cf^ and 
compared them with the relative variation in the C;'s induced when changing 
ttHe and ^He by a 50%. In all of the cases the first quantity is several orders 
of magnitude greater than the variation of the C^'s. Therefore, in order to 
analyze WMAP5 data, there is no need to modify these parameters. 

To put constraints on the variation of a and rrie during recombination 
time in the detailed recombination scenario studied here, we introduced the 
dependencies on the fundamental constants explicitly in the latest version 
of Recfast code [18], which solves the recombination equations. We per- 
formed our statistical analysis by exploring the parameter space with Monte 
Carlo Markov chains generated with the publicly available CosmoMC code 
of ref. [19] which uses the Boltzmann code CAMB [20] and Recfast to 
compute the CMB power spectra. We modified them in order to include 
the possible variation of a and rrie at recombination. We ran eight Markov 
chains and followed the convergence criterion of ref. [21] to stop them when 
R — 1 < 0.0180. Results are shown in Tabled] and Figure [21 

The observational set used for the analysis was data from the WMAP 
5- year temperature and temperature-polarization power spectrum [6], and 
other CMB experiments such as CBI [22], ACBAR [23], and BOOMERANG 
[2H [25] . and the power spectrum of the 2dFGRS [26]. We have considered 
a spatially-flat cosmological model with adiabatic density fluctuations, and 
the following parameters: 



where flcDAih'^ is the dark matter density in units of the critical density, 9 
gives the ratio of the comoving sound horizon at decoupling to the angular 
diameter distance to the surface of last scattering, r is the reionization optical 
depth, Hs the scalar spectral index and is the amplitude of the density 
fluctuations. 

In Fig. [2] we show the marginalized posterior distributions for the cosmo- 
logical parameters, Aa/a^, and Ame/(me)o, which are the variation in the 
values of those fundamental constants between recombination epoch and the 
present. The three succesively larger two dimensional contours in each panel 
correspond to the 68%-, 95%-, and 99%- probability levels, respectively. In 
the diagonal, the one dimensional likelihoods show the posterior distribution 
of the parameters. 



P= [QBh^ncDMh\e,T, 



Aa AvUf. 




(13) 
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Figure 2: Marginalized posterior distributions obtained with CMB data, in- 
cluding the WMAP 5-year data release plus 2dFGRS power spectrum. The 
diagonal shows the posterior distributions for individual parameters, the 
other panels shows the 2D contours for pairs of parameters, marginalizing 
over the others. 

In Table [1] we show the results of our statistical analysis, and compare 
them with the ones we have presented in ref. [4|, which were obtained in the 
standard recombination scenario (i.e. the one described in [27], which we 
denote PS), and using WMAP3 [28l [29] data. The constraints are tighter 
in the current analysis, which is an expectable fact since we are working 
with more accurate data from WMAP. The bounds obtained are consistent 
with null variation, for both a and mg, but in the present analysis, the 68% 
confidence limits on the variation of both constants have changed. In the 
case of a, the present limit is more consistent with null variation than the 
previous one, while in the case of mg the single parameters limits have moved 
toward lower values. To study the origin of this difference, we perform an- 
other statistical analysis, namely the analysis of the standard recombination 
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Table 1: Mean values and la errors for the parameters including a and me 
variations. NS stands for the new recombination scenario, and PS stands for 
the previous one. 

scenario (PS) together with WMAP5 data and the other CMB data sets and 
the 2dFGRS power spectrum. The results are also shown in Table [1] We 
see that the changes in the results are due to the new WMAP data set, and 
not to the new recombination scenario. In Fig. [3] we compare the probability 
distribution for Aa/a^ and also for Ame/(me)o, in different scenarios and 
with different data sets. 

In Fig. Hlwe compare the 95%- probability contour level for the parame- 
ters, and their one dimensional distributions, for two different analysis in the 
standard recombination scenario, namely the one with WMAP5 data (dashed 
lines) and the one with WMAP3 data (solid lines). The contours are smaller 
in the former case, which is expectable since that data set is more accurate. 
For the fundamental constants, the contours notably shrink. Moreover, the 
constraints are shifted to a region of the parameter space closer to that of null 
variation in the case of a. On the other hand, limits on the variation of nif. 
are shifted to negative values, but still consistent with null variation. From 
the one dimensional likelihoods we see that the the peak of the likelihood has 
moved for VLbh?. The obtained results for the cosmological parameters are 
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Figure 3: One dimensional likelihood for ^ (upper row) and (lower 
row). Left: for WMAP5 data and two different recombination scenarios. 
Right: comparison for the standard recombination scenario, between the 
WMAP3 and WMAP5 data sets. 

in agreement within Icr with the ones obtained by the WMAP collaboration 
|30j . without considering variation of fundamental constants. 
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